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Abstract. One of the nicest results in cosmological perturbation theory is the 
analytical resummaton of the leading corrections at large momentum, which was 
obtained by Crocce and Scoccimarro for the propagator in [1]. Using an exact 
evolution equation, we generalize this result, by showing that a class of next-to- 
leading corrections can also be resummed at all orders in perturbation theory. The 
new corrections modify the propagator by a few percent in the Baryonic Acoustic 
Oscillation range of scales, and therefore cannot be neglected in resummation schemes 
aiming at an accuracy compatible with future generation galaxy surveys. Similar tools 
can be employed to derive improved approximations for the Power Spectrum. 
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1. Introduction 

Cosmological perturbation theory (hereafter, PT, for a review, see [2]) has attracted a 
renewed interest in the recent literature. On the one hand, future galaxy surveys will 
measure - at the percent level- the statistical properties of matter distribution in a 
range of scales and redshifts in which linear perturbation theory (LT) predictions for 
the power spectrum (PS) and higher order correlation functions are not more accurate 
than 0(10%). On the other hand, due to practical limitations in machine-time, N- 
body simulations are far from being the ideal tool to scan over parameters and models, 
and involve larger complexity when extra components (such as massive neutrinos, 
a non-minimally coupled quintessence field, ...) or non-gaussian initial conditions 
are considered. Moreover, a series of works, pioneered by Ref. [I] by Crocce and 
Scoccimarro (hereafter, CS), has successfully investigated the possibility of improving 
PT by "resumming" perturbative contributions at all orders. Different resummation 
schemes have been proposed and their results for the PS in the ACDM cosmology have 
been compared to N-body simulations in the Baryon Acoustic Oscillations range of 
scales [SlIllElElEllHlinilinillllinillSlElIIS], showing in most cases an agreement in 
the few percent range. These methods have been applied also to less standard -although 
interesting - cosmological scenarios, such as mixed dark matter (i.e. cold dark matter 
+ massive neutrinos) [IB] , coupled dark matter-dark energy models [TF], and models 
with various types of primordial non-Gaussianity [T8] . 

The main reason for the success of these resummation methods, compared to 
traditional PT, was nicely discussed in CS. In PT, each successive order in the expansion 
in powers of the linear PS becomes more and more important at small scales, or large 
Fourier momenta, k. Therefore, the perturbative expansion unavoidably breaks down 
at high fc's and low redshifts. On the other hand, the leading contributions at large fc's 
have a simple form which, in some cases, allows them to be summed to all orders, giving 
a convergent result. This was shown analytically for the propagator, i.e. the - properly 
normalized - cross-correlator between the final density, or velocity, perturbations and 



the initial ones (for its exact meaning, see Eq. (26)). In this case, the large- momentum 
leading order contribution, at n-th order in PT, grows as (— l) n k 2n , but the sum exhibits 



a nice Gaussian damping, cx exp(— k 2 a 2 D 2 /2), with D the linear growing mode, and o 2 



the one-dimensional dispersion velocity (see Eq. (36)). 

In other words, the bad UV behavior of PT can be greatly alleviated if one 
reorganizes the series expansion, for instance using as a zero order approximation the 
resummed propagator of CS instead of its linear approximation. This idea forms the 
basis of the 'renormalized perturbation theory' (RPT) approach [19]. Similar analytical 
results hold also for a special class of higher order correlation functions, the 'multi-point 
propagators considered in [201 [21]. On the other hand, for the most directly observable 
quantities such as the PS, the bispectrum, and higher order correlators, the resummation 
program cannot be carried out analytically, and the semi-analytical methods are needed 

^EiiEicziEiiiniiniiiaiiainiiE!. 
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In this paper, we will discuss how to go beyond the CS result for the propagator 
by taking into account subleading corrections which were neglected in the analytic 
resummation described above. We will identify a class of next-to-leading contributions 
and we will show that the resummation performed in CS can be generalized to include 
this larger class of perturbative corrections. Our results indicate that the new resummed 
propagator corrects the CS one by a few percent in the BAO range, and therefore the 
new contributions should be taken into account. 

The basic tool to go beyond the CS resummation is provided by exact evolution 
equations for the propagator, which we will derive. We will discuss in which 
approximation the solution of this equation gives the CS propagator and then explore 
two different schemes to incorporate next-to leading corrections. 

This paper is organized as follows. In Section [2] we review the derivation of 
PT in the useful compact form introduced in [19j. In Section [3] we introduce the 
generating functional for the statistical correlators, along the lines of Ref. [6], recall 
the diagrammatic language useful to discuss PT and its resummations, and derive the 
exact evolution equation for the propagator. In Section[4]we clarify the relation between 
the evolution equation and the resummation of 'chain-diagrams' performed in CS, and 
then, in Section [5] we discuss how to go beyond the CS result by taking into account the 
renormalized PS in the chain-diagrams. In Section [6] we present our numerical results 
and, finally, in Section [7] we discuss them and give our conclusions. 



2. Nonlinear fluid equations and the Propagator 

We can write the three fluid equations (continuity, Euler, Poisson) in Einstein-de Sitter 
cosmology as follows 

^ + v-[(i + <Uv] = o, 

or 

V Hv + (v • V)v = -V0, 



Or 



3 



2 

where "H = d log a/dr is the Hubble Parameter in conformal time, while £ m (x, r) and 
v(x, r) are the DM number- density fluctuation and the DM peculiar velocity field, 
respectively. 

Defining, as usual, the velocity divergence 6*(x, r) = V • v(x, r), and going to 
Fourier space, the equations in ([!]) give 

^ k ^ + 0(k,r) 
or 

+ J d 3 q d 3 p 5 D (k-q- p)a(q, p)0(q, r)5 m (p, r) = , 

afl M + ^(k 1 r) + ^l(k,r) 
or 2 

+ Jd 3 qd 3 p6 D (k-q-p)f3(q,p)9(q,T)9(p,T)=0. (2) 
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The nonlinearity and non-locality of the fluid equation are encoded in the two functions 

/ n (p + q)-q at n (p + q) 2 p-q ^ 

which couple different modes of density and velocity fluctuations. 

One can write Eqs. Q in a compact form pp. First, we introduce the doublet tp a 
(a = 1,2), given by 

( <PlQz,V) \ = -r, ( MM) \ m 

where the time variable has been replaced by the logarithm of the scale factor, 

i a 

rj = log — , 

ai n being the scale factor at a conveniently remote epoch, in which all the relevant scales 
are well inside the linear regime. 

Then, we define a vertex function, 7 a b c (k, p,q) (a, b, c, = 1,2) whose only 
independent, non- vanishing, elements are 

7i2i (k, p, q) = - 5 D (k + p + q) a(p, q) , 

722 2 (k, p, q) = 5 D (k + p + q)/3(p,q), (5) 

and 7i2i (k, p, q) = 7n 2 (k, q, p). 

The two equations ^ can now be rewritten in a compact form as 

d v v? a (k, 77) = -Q a b Vb{K 77) + e r? 7 a6c (k, -p, -q)(p b (p, 77) <p c (q, 77), (6) 

where 

( 1 ~M 

Q = (7) 
3 3 V ; 

V "2 2 / 

Repeated indices are summed over, and integration over momenta q and p is understood. 

To extend the validity of this approach to ACDM, we will reinterpret the variable 
rj as the logarithm of the linear growth factor of the growing mode, i.e. [21 [fl [9], 

77 = HD/D m ) , (8) 

and we redefine the field in Eq. Q as 

/ <pi(k,rj) \ = ( 5 m (k,77) \ 

{^(Krj))- {-e(Kv)/nf) ' lJ 

with / = d In D/d In a. As discussed in j[9], the above approximation is accurate at 
better than 1% level in the whole range of redshifts and scales we are interested in. 

If we consider the linear equations (obtained in the e^^iabc limit) we can define 
the linear retarded propagator as the operator giving the evolution of the field (p a from 

Vin to 7/, 

fa ( k > V) = 9ab{V, Vin)fb(K Vin) ■ (10) 
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The linear propagator obeys the equation 

(S a bd v + Qab)gbc(V, Vin) = ^acM 7 ? - 7] in ) . (11) 

with causal boundary conditions. It is given explicitly by the following expression pQ , 
9ab(Va,Vb) = [B + Ae~ 5 W" b )] 0( _ Vb) (12 ) 

L J ab 

with 9 the step-function, and 

B=if 3 2 ) and A=i( 2 " 2 ) . (13) 

5 v 3 2 ; 5 v -3 3 y y J 

The growing (<^ a oc const.) and the decaying (ip a oc exp(— 5/2rj a )) modes can be selected 
by considering initial fields (p a proportional to 

Ua= {\) Va = ( -3/2 ) ' (14) 

respectively. 

3. The non-linear propagator and its evolution 

Following the path-integral formulation of cosmological perturbation theory introduced 
in [H] we can derive exact evolution equations for the propagator and for the PS. In this 
section we will give a short review of the formalism introduced and discussed in [6] , and 
we will obtain the exact evolution equation for the propagator, which will be solved in 
various approximations in the remaining sections. 

The generating functional for the time-dependent correlators between perturbations 
is given by 

Z[J a , K b ; P°] = J T>(p a Vxb expj-i J dr] a dr] b x a Pab S (Va)S(r] b )xb 

+ i J drj [xa9db<Pb - e v J a bcXa<Pb<pc + J a <Pa + KbXb] J , (15) 

where J a and K b are sources for cp a and Xb respectively, and Pab(k) is the PS at the 
initial time r\i n = 0. In deriving the above expression we have assumed Gaussian initial 
conditions. Non-Gaussian initial conditions can be taken into account by including a 



non- vanishing bispectrum, trispectrum, etc., in the first line of Eq. (15). Derivatives of 



Eq. (15) w.r.t. the sources J a and K b give all the possible statistical correlators involving 
the fields ip a and Xb- As usual, it is more convenient to discuss connected correlators, 
which can be derived from the generating functional, 

W = -i\ogZ. (16) 
In the following, we will be interested in the PS 

{ip a (K i] a )ip b (k', rjb)) = 8 D (k + k')P ab (k; rj a , rj b ) , (17) 
and in the propagator, 

(ip a (k,T] a )xb(k',Vb)) = iS D (k + k')G ab (k;ri a ,r] b ) . (18) 
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They are given by the second derivatives of W, according to the relations, 
5 2 W 



5J a 5J b 
S 2 W 



6J a SK b 
8 2 W 



SK a 5J b 
S 2 W 



Ja,K b =0 



Ja,K b =0 



Ja,K b =0 



i5 D (k + k')P ab . 



-<5 D (k + k')G ab , 
-5 D (k + k')G ba , 



SK a 6K b 



(19) 



Ja,K b =0 



One can also consider the effective action, built in terms of the average fields 



-rr KMX, ,- bW[J c ,K d ] SW[J c ,K d ] 

(fa [Jc, K d (k, rf) = — , Xa [J c ,K d ](k, 77) = — , (20) 

OJa 0K b 

where the functional derivatives are evaluated at generic {i.e. non-vanishing) values for 
the sources J and K. The effective action is given by the Legendre transform of W, 

F[<p a , Xb] = W[J a , K b ] - Jd v d 3 k {J a (p a + K bXb ) , (21) 

and its derivatives with respect to <p a and Xa give rise to the one-particle irreducible 
Green functions (1PI). The two-point 1P1 functions are given by 

s 2 W a ,x b ] 



Sip a 5ip b 

s 2 r[^ a , Xb } 



0. 



S<p a Sxb 

s 2 r[^ a , Xb } 



SXaSXh 



<Pa,Xb=0 



<Pa,Xb=0 



<Pa,Xb=0 



<Pa,Xb=0 



gj -S a 6j 5 D (k a + k b 



(dba -^ba)S D (k a + k b ) , 

iP° b (k)6(r])5(7] b ) + i$ ab ) 5 D (k a + k b 



(22) 



where we have isolated the 'free' parts, g ab l and P® b , from the 'interacting' ones, S ab and 

The 1PI two-point functions in Eq. (22) are related to connected ones in Eq. (19) 
by the following relations, 



P a b(k; r) a , r] b ) = G ac (k; r) a , 0)G bd (k; rj b , 0)P° d (k) 

fVa pb 

+ dsi ds 2 G ac (k;r] a ,s 1 )G bd (k;ri b , s 2 )$ c d(k; s 1: s 2 ) , 
Jo Jo 

(23) 



and 



-1 



Gab(k;Va,Vb) = 9 ba - S fea (k;Va,Vb) 



(24) 
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where the last expression has to be interpreted in a formal sense, that is, 

G a b(k; Va,Vb) = gab(Va,Vb) 



+ 



'la 



'lb 



dsi \ ds 2 g ac (rj a ,s 1 )T, cd (k; s 1 ,s 2 )g db (s 2 ,rj b ) + 



(25) 



While the physical meaning of the PS - namely the cross-correlator of the density- 
velocity fields ip a and tp b computed at times r\ a and r\ h - is clear from its definition, 
Eq. (17), the presence of the auxiliary field x m the definition (18) makes the 
interpretation of the propagator more obscure. Insight on its physical meaning can 
be obtained from Eq. (23), by sending one of the two times, e.g. rjb to the initial time, 
i.e. rjb ->■ + , 

P ab (k; Va , 0) = G ac (k; Va , 0) + f" ds G ac (k; Va , s) $ cb (k; s, 0) , (26) 

Jo 

where we have used the property of the full propagator G b d(k;rj,0) — > 6 b d for r\ — > + , 
see Eqs. (12) and (25). The second term at the RHS of the equation above vanishes 
if the statistics is Gaussian at the initial time, therefore, in this case, the propagator 
connects the initial PS, Pa b (k) = P ab {k;rj a = rjb = 0), to the cross-correlator between 
the initial fields and the 'final' ones evaluated at rj a > 0. Taking the initial PS on the 
linear growing mode, 

P° ab (k) ~ P°(k)u a u b , (27) 



where the u a vector has been defined in Eq. (14), we see that the late-time cross 
correlators are entirely given in terms of P°(k) and the two combinations, 

G a (k;r!,0)=G ac (k;r!,0)u c , (a = 1,2), (28) 

which will be considered in the following. 

We notice here that Eq. (26) holds in general, i.e., even when the initial time is 
taken to correspond to a lower redshift, where non-linearities and non-gaussianities have 
to be taken into account. We will make use of this relation in sect. [U 

In general, all kind of N-point correlators, and in particular the two-point functions, 
can be computed in perturbation theory (PT). From the path integral formulation 
discussed above one can derive Feynman rules for the building blocks of PT, namely, 
the free propagator, the linear PS, and the interaction vertex [5], from which all the 
higher order correlators can be built. These are summarized in Fig. [T] continuous and 
dashed lines indicate (p a and Xa fields, respectively. 

Equations (23) and (24) by themselves, however, do not rely on PT, and therefore 
offer the opportunity of computing the two-point correlators non-perturbatively. A 
convenient way to deal with the propagator in a non-perturbative way is to cast Eq. (25) 
in a closed form [19] [f] 

G ab (k; rj a ,rj b ) = g ab {r, a ,r, b ) 

rr\a rs\ 

+ / ds x \ ds 2 g ac (Va, Si)E cd (fc; s l} s 2 ) G db (k; s 2 ,r) b ) , 

Jr\b J Tit, 



| In turbulence theory this result is well known, see for instance Ref. 
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propagator: -ig ab {rj a^r\h) 



a power spectrum: P° b (k, r? a ,??b) 



i 

A 



interaction vertex: —ie v j a bc(k a , k fc , k c 



Figure 1. The Feynman Rules for cosmological perturbation theory 



(29) 



which is equivalent to (25), as can be shown by expanding the full propagator G at the 



RHS iteratively in the 'self-energy' E. Then, by deriving Eq. (29) with respect to r] a we 
get 

d Va G ab (k; r) a ,r}b) = -Q ac G cb (k; r] a ,r] b ) 

+ [ dsS ac (/c; 7] a ,s)G cb (k; s,rj b ) , (30) 
J rib 



which gives the exact {i.e. non-perturbative) evolution equation for the full propagator, 
and will be the starting point for our evaluation of G ab . 



4. Factorization and the Crocce-Scoccimarro propagator 



In order to get insight on the content of the exact evolution equation, Eq. (30), we will 
consider a perturbative expansion for the propagator and the 'self-energy', 



G ab {k; Va,Vb) = G ri( k > Va,Vb) , 



71=0 

OO 



(31) 



£»&(&; Va,Vb) = H S afe Va,Vb), 
n=l 

where, as usual, the index n counts the number of power spectra contained in the n- 
th order contributions to G ab and S ab . Notice that at zeroth order E ab receives no 
contribution, while 



G { a b {k; r] a , rj b ) = g ab {rj a - rj b ) . 



(32) 



Inserting (31) in (30), and equating terms of the same order, we get the evolution 



equation for the n-th order contribution to the full propagator, 

d Va G { a b \k\ Va,Vb) = - ^l ac G^ { k ] T] a ,T] b ) 
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Figure 2. sic , the 1-loop contribution to the self-energy. 




Figure 3. Chain-diagrams. 



n-l 



'la 



+ e n , £ / ds E&-'> (A;; 7; a ,s) G^ b '(k; s, rj b ) , (33) 
i=o ^ 

where nj o is zero for n = and one otherwise. 

In the large external momentum limit the leading diagrams contributing to the last 
line of Eq. (33) are the so-called chain- diagrams (see Fig. [3]) already discussed by CS 
in [1], all the other contributions being suppressed by inverse powers of k. As discussed 
in detail in Appendix A in this limit the sum can be computed analytically, giving the 
remarkable factorized result 



£ rdsEf j \k;r ]a ,s)G ( £(k; s, Vb 

G^ik; Va , Vb ) rds^}{k 



large k ^-y(n-l) 



(34) 



where we used the boldface for the index a in the last line to indicate that it is not 
summed over. represents the 1-loop contribution to the self-energy corresponding 
to the diagram Fig. [2} In the large momentum limit one gets 

ds sW(fc; r) a , s) u c la -^> k -JfeV e" a (e^ - e % ) , (for a = 1, 2) . (35) 



Va 
'lb 



with 



a' 



P°(q) 



Q 1 



(36) 



Summing Eq. (33) over n, we get the evolution equation for the full propagator in 



the large momentum limit 

d Va G ab (k; r] a ,rj b ) = -fi ac G c6 (/c; r} a ,r) b ) 

+ G ab (k; i] a , r) b ) i ds E^(k; i] a , s) 

Jrih 



(37) 
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At large k the differential Eq. (37) can be easily integrate in rj a and it yields 

G ab (k; 7] a , r) b )u b = exp -k 2 a 2 (for a = 1, 2) , (3£ 



this reproduces the Gaussian decay of the large momentum propagator found by CS in 
Ref. pQ 

In the opposite limit, k 0, higher order contributions to the propagator are 
suppressed, and linear perturbation theory is recovered. In order to take into account 
the first non- linear corrections in this limit, one can truncate the series in (31) at n = 1, 



i.e. at 1-loop order. Moreover, we will consider the evolution of the two combinations 



G a , introduced in Eq. (|28|). Therefore, the relevant term in the sum of Eq. (33) in the 

(39) 



small k limit is the one for n = 1, namely, 

fVa , . pa , . 

/ dsE ( J(k;r]a,s)g cb (s-r] b )u b = dsE ( J(k; r) a ,s)u c , 

•'Vb JVb 

where we have used the property of the linear propagator, 

g a b(v)u b = u a . 

Modulo terms at least of 2-loop order, the above expression can be rewritten as 
G ab (k; rj a , rj b ) u b / ds T, { J(k; rj a ,s)u c , 



(40) 



(41) 



which gives the same equation as the one for large k, Eq. (37), contracted by u b § It is 



remarkable that the same factorization holds in the two limits of large and small k. 
For any k value the 1-loop contribution to the self-energy reads 

Va: s) = 

Ae Va+s J d 3 q'j ade (k, -q, q - k)P°(q)u d u f g eh (r] a - s)7 fe/c (k - q, q, -k) , 

(42) 

and allows us to compute the propagator both in the small and large momentum limits 



with the evolution Eq. (37). Notice that the PS appearing in Eq. (42) is the linear one 



The above discussion clarifies the comparison between the present approach and the 



one in [TJ. Indeed, both give the same result, Eq. (38) for k — > oo, and both reproduce 
the 1-loop propagator for k — > 0. For intermediate k's the two approaches give different 
prescriptions to interpolate between their common limits. In the present section, the 



interpolation is based on two approximations. First, we factorize the integral of Eq. (30) 

as 

ds S ac (/c; r] a , s) G cb (k; s, r] b ) ~ G ab (k; r) a , r) b ) / ds £ ac (£;; r] a ,s)u c , 

Vb Jr\b 

(43) 

§ If we do not contract by Ub, the factorization still holds exactly for the individual components of 
the propagator at small k if one takes the 775 — > —00 limit. For finite 775 , the factorization is not exact 
anymore, but it is anyway a very good approximation. 



Next-to-leading resummations in cosmological perturbation theory 



11 



where S ac refers to the approximation chosen for the two-point 1PI function. Second, 
we consider the purely 1-loop 'self-energy' in the factorized expression at the RHS above 



SacO; 77a, S) - ^ (h, 7] a , s) 



(44) 



These are the two approximations needed to pass from the exact equation (30) to the 
approximated one, Eq. (37), which we will solve for any k. On the other hand, in 
PQ the interpolation is obtained by advocating an exponentiation procedure. We can 
directly check the differences between the two procedures by taking the rj derivative of 
the propagator given in Eq. (41) of [TJ. In the limit r)i — > — oo, the comparison simplifies 
considerably, and one finds an evolution equation given by Eq. (37) plus an extra term 
at the RHS, given by 



where 



U ab 



exp (G^JucJ - G^Uc exp (G^u d 



f ds f ds'g ac (r] - s)T,^(k; s, s')g db (s' 



(45) 



(46) 



is the 1-loop contribution to the propagator, and il a b has been defined in Eq. (|7|). The 
extra term is at most of two-loop order and it vanishes both in the large and small 
momentum limits, as it should. Therefore, in the rest of this paper, we will refer to the 
propagator obtained by solving Eq. (30) in the double approximation (43) and (44) as 
the CS propagator, even though it differs from the one of ref. p] by subleading terms at 
intermediate k. We stress again that the difference between the propagator computed 
along the lines described in this section and the one in [TJ is entirely due to the different 
interpolation procedures, while the classes of diagrams resummed are exactly the same. 

In the following we will go beyond the approximation of this section in a consistent 
way. As mentioned above, we will keep the factorized form of the equation for any k, 
but will consider new contributions to the 'self-energy', therefore improving over the 
approximation of Eq. (44). As we will show, this corresponds to resumming a larger 
class of diagrams than just the infinite chains considered by CS. 



5. Extended factorization: the renormalized chain-diagrams 



The large-A; factorization property, Eq. (34), holds for a more general class of diagrams 



than the 'chain'-ones discussed by CS. This is the main result of this paper, and is 



proved in Appendix A There, we show that by replacing all the linear power spectra, 
(q) = P°(q)u a Ub, appearing in the chain-diagrams by a - for now - generic non-linear 



PS, of the form Pa b l (q; s a , s b ), one still gets Eq. (34) in the large-/c limit, where now, 

t ac (k;Va,s)^Zl^(k; Va ,s). (47) 

S^ nl is obtained from the 1-loop self-energy by replacing P® b by P£ b l , as indicated in the 
diagram in Fig. |4j and has the large-momentum limit 



Va 



dsZ p J»\k- Va ,s)u t 



large k 



'lb 




ds e"°" 



d 3 q 



. (4* 
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The effect of the inclusion of these subleading corrections is clear: at large momenta the 
propagator still decays exponentially, but with the decay law of Eq. (p8| replaced by 



G ab (k; r] a ,r] b )u b = exp —k a nl (r] ai r] b ) 



una - e vi)r 



(for a = 1,2), (49) 



where 



'lb 



*1 



dsi / ds2 e 



S1+S2 



lib 



:| ^22(9; s l> s 2) 



d*q 



Q 



(50) 



Notice that only the "22" (i.e. velocity-velocity) component of the PS appears in the 
exponential above. Since it is known that this component receives negative corrections at 
the non-linear level, (see, for instance, [9j[13]), we expect that the improved propagator 
will be enhanced w.r.t the CS one at large k. 

In the opposite limit, k — > 0, we are no longer guaranteed that Eq. (47) is still a 
good approximation, and that the contributions to the S ac obtained by replacing the 
linear PS with the non-linear one are the only leading ones. In other words we consider 
again the n-order contribution expressed in Eq. ( 33 ) where now the upper indices count 
the number of non-linear PS in a given contribution to X and G. In this context, in 
the small-fc limit, again the relevant term in the sum of Eq. (33) is the one for n — 1. 
Following the same argument of the previous section the non linear part of Eq. (33) 
becomes 



G ab {k; rja, Vb) u b [ ds T^ nl (k; n a , s) u c 



(51) 



getting the factorization property in the small- limit. At this level we avoid any 
double counting problem, however the application of this perturbative criterion must be 
analyzed carefully by comparing the results with the standard perturbative computation. 
In this respect, a powerful guiding criterium is the requirement that linear theory is 
recovered for small momenta. Indeed, from the exact equations (25) and (29), one 
concludes that, in order to have G ab (k; r] a ,7]h) — > g a b(Va — Vb) as k — > the 'self-energy' 
has to vanish in this limit 

s a #;^ s)^o. (52) 

therefore we should also have £^ c Snl (A;; r] a , s) vanishing for k — > 0. In the approximation 
(51) this is not automatically realized. Indeed, one finds 

rva 



r a dsY? a ^\k- Va ,s) 

J Vb 



Uc 



8 a i r ds eVa+S I d3( l \92d(Va - s)P^(q; r] a , s) - g ld (r] a - s)P^(q; rj a , s) 



Vb 



(53) 

Of course, if one puts back the linear PS, P® b = P°u a u b , in place of P£ b l in the 
expression above, one recovers the 1-loop self-energy, which vanishes in the A; — ?• 
limit as can be directly checked from Eq. (53), using Eq. (40) (the first non- vanishing 
contribution goes as k 2 ). 
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Figure 4. S^f nl (A:; rj a , s): 1-loop self-energy with the insertion of the non linear PS. 

... + „-^._^.. + . 

[~| = non linear PS 

Figure 5. Renormalized chain-diagrams. 

On the other hand, moving a step further and using the 1-loop result for P^, 
namely, including the diagrams of Fig. [6] in the computation of the self-energy, one finds 
a non- vanishing limit for k — > 0. Indeed, one can check that the contributions from the 
remaining 2-loop diagrams listed in Fig. [7] exactly cancel those of Fig. |6j recovering in 



this way the physical requirement of Eq. ( 52 ) . 



At large k the contributions V and VI in Fig. [7] would give rise to chain-diagrams 
for the propagator with the insertion of a linear PS. These contributions are already 
taken into account by diagram I in Fig. |6j so, in order to avoid double counting, we do 
not have to include them. The remaining diagrams, VII-IX, are subdominant at large k 
w.r.t. the chain-diagrams. On the other hand, at small k, all the diagrams in Fig. [7] are 
essential in order to recover linear theory. Therefore, an approximation to S ac giving 
the '1-loop renormalized' chain-diagrams (i.e. the chain-diagrams with the 1-loop PS 
replacing the linear one) in the large k limit, and recovering linear theory for k — > 0, is 
given by 

t ac (k; r]a,s)- K! U (k; r? a , s) + hmE^(A;; r] a , s) 

fc-»0 

= K! U (k; Va ,s)-KmZ™(k; Va ,s), (54) 
fc— >o 

where S^ 11 and S^ c rcst are the contributions to the self-energy computed from the 
diagrams of figs. [6] and [TJ respectively. The S^ c rcst (A;; r\ a , s) term gives also contributions 
of order k 2 when k approaches to zero, that we do not include with the prescription in 



(54). Therefore, our resummed propagator reproduces the linear one in the k — > limit, 
with the 0(k 2 ) terms given by the 1-loop one plus some, but not all, of the 2-loop and 
higher contributions. 



Finally we point out that using the 1-loop PS corresponds to compute Eq. (33) for 



n = 1 and for n = 2 with the usual non-renormalized Feynman diagrams. Therefore 
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Figure 6. T,^ n (k; r\ a ,s): the self-energy with the insertion of the PS up to 1-loop 
order. 
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Figure 7. T,^ BBt : the 2-loop contributions to the self-energy not included in E^f 11 . 



one can check that the approximation (51) holds in this case. Indeed computing the 



non linear term of Eq. (33) at 2-loop level one gets 
rva 



ds Va , s) G ( ^(k; s,rj b ) + rj a , s) g cb (s - n b ))u b . (55) 



'ib 



In Eq. (55) the first term goes as k A when k approaches to zero and results subdominant 
w.r.t. the second one (that goes as k 2 ). This allows to consider just the second term and 
to advocate the procedure involved in the previous section to achieve the factorization 
property also at 2-loop order. This argument justifies our factorization procedure given 



by expression (51) also for small values of k. 



A further extension of the resummation program is to use as non-linear PS in the 
computation of S^f nl the one computed with the Time Renormalization Group (TRG) 



approach introduced in [9] and briefly reviewed in Appendix B As discussed in [9], the 
TRG equations, truncated at the bispectrum level, incorporate perturbative corrections 
in which the PS lines are iteratively replaced by their 1-loop corrections. This procedure 
resums perturbative contributions at all orders, of which some 3 and 4-loop examples 
are listed in Fig. [8j The TRG approach is able to reproduce the non-linear PS at z = 
up to k ~ 0.2 h/Mpc (that is, in the BAO region) at the few percent level [HI |2^] . 
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+ 



Figure 8. Some of the contributions included in E a b when the non-linear PS is given 
by the one computed via the TRG. 



The TRG, as it is, gives the non-linear PS computed at equal times, that is 
(g;?7, 77). On the other hand, in the computation of the self-energy, we need the 



PS computed at different times, see Eqs. (48, 53). The relation between the equal-times 



PS and the one computed at different times can be read from Eq. (26) where we now 
take the initial time rj = to be any generic time, 77 = s. If s corresponds to a low 
redshift, non-linearities and non-gaussianities cannot be neglected. They are encoded 
in the non-linear initial PS, that we will compute with the TRG, and in the irreducible 
function $ c d(/c; rj a , s), which gives a subdominant contribution that we will compute in 
PT. When dealing with Pj b RG (q; r] a , s) we will therefore replace it by 



G ( ^ c s (q ] r ]a ,s)P] b KG (q-s,s) 

+ J e Sl J d 3 qg ac (r) a - Si)7 cde (k, -q, q - k) 
x 9df(s 1 - s) g eg {si - s) B bfg (k, -q, q - k; s) 



(56) 

where, for practical reasons, we have replaced the full propagator with, G^ c , the 
propagator computed a la Crocce-Scoccimarro according to the approximation discussed 
in Section [4], and B is the tree-level Bispectrum computed via perturbation theory. 

As in the computation of the 11 , we find a non- vanishing k — > limit for the self- 
energy computed with p TRG . We follow the same arguments discussed above and, also 
in this case, we incorporate the relevant corrections in this limit by using the prescription 



Z ac (k; Va,s) ^ Sf™ G (/c; T] a , s) - lim E£~*(fc; r] a , s 



pTRG 



(57) 



To summarize the results of this section, the evolution equation for the 
propagator including next-to-leading corrections (the renormalized chain-diagrams plus 
the contributions needed to recover the proper k — > limit) is given by 

dr, a G ab (k; r] a ,r] b ) = - Vt ac G cb {k- } r] a ,r} b ) 

+ G ab (k] Va,Vb) ds£ac(fc; Va,s)u c , (58) 



</{, 



which we will solve for any scale k using the two approximations, Eqs. (54) and (57) 
for the self-energy E ac . 




Figure 9. The density (left panel) and velocity (right panel) propagators at redshifts 
z = and z = 1. The black-dotted lines are the propagator computed in the CS 
approximation. The purple-dashed lines are obtained using the 1-loop PS in E, while 
the blue-continuous lines are obtained by using the PS from the TRG. 



6. Results 

We investigate a ACDM cosmology close to the best-fit model (fl m = 0.25, Qbh 2 = 
0.0224, h = 0.72, n = 0.97 and a 8 = 0.8). The initial time, rj = 0, is taken to 
correspond to the physical redshift Zi n = 100. At 77 = we set the initial conditions for 



the evolution equation (58) and for the TRG equations needed to compute the PS (see 



Appendix B). We set the initial conditions for the PS by matching it with the linear PS 



obtained by the CAMB code [25J. For the propagator, we use the explicit expression 



given in Eq. (12). The integration of the TRG equations requires also initial values for 
the bispectra. We set them to zero, i.e., we neglect all non-Gaussianities generated at 
redshifts higher than z = 100. 



In Fig. [9] we plot our results for the propagators G a (k; r) a , 0) defined in Eq. (28) (G\ 
in the left panel, G2 in the right panel), computed at final times rj(z) corresponding to 
redshift z = and z — 1. The dotted black lines correspond to the Crocce-Scoccimarro 



result, i.e. to the integration of Eq. (37), where the linear PS has been used to compute 
E^. The dashed red lines are obtained by using the 1-loop approximation for the 
non-linear PS in E ac , while the continuous blue lines are obtained by using the TRG 
PS. 



In Fig. 10 we plot the relative difference between the propagators computed with 
the two different approximations for the non-linear PS and the one computed in the 
Crocce-Scoccimarro approximation. 

As a general trend, the effect of the inclusion of the new class of diagrams considered 
in this paper leads to a weaker damping of the propagators at intermediate and large 
fc's, compared to the one obtained considering only the chain-diagrams of CS. The effect 
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Figure 10. Relative differences between the improved propagators and the one 
obtained in the CS approximation. Line-codes as in Fig. [9] 
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Figure 11. The dependence on the UV-cutoff of the improved propagator obtained 
by using the 1-loop approximation for the non-linear PS. Dashed lines are obtained 
using q uv — 1 h/Mpc -1 , while dash-dotted lines are obtained with q uv — 2 h/Mpc -1 . 



is stronger for the velocity propagator than for the density one. 

The computation with the 1-loop PS suffers from an intrinsic uncertainty, due to 
the dependence on the (UV) momentum cutoff employed in the loop integrals. This 
is an unavoidabe limitation of the 1-loop approximation, due to the fact that the 1- 
loop PS takes unphysical negative values at large g's, especially when the two time 
arguments are very different. In Fig. 11 , we show the cutoff dependence by plotting the 
same quantities as in Fig. 10, computed using the 1-loop approximation for the PS and 
UV cutoffs g max = 1 and 2h/Mpc. The cutoff dependence is quite strong for the G2 
propagator, showing that its computation using the 1-loop PS is clearly unreliable at 
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Figure 12. ajff a at z = 1 computed from G a (see Eq. 59). Line-codes as in Fig. &9 
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Figure 13. The same as Fig. 12 but at z = 



low redshift. On the other hand, the results for Gi j2 obtained using the PS computed 
with the TRG, which is always positive, do not exhibit UV problems. In this paper 
we use always the g max = 1 h/Mpc cutoff limiting in this way the effects due to the 
unphysical behavior of the 1-loop power spectrum approximation. 

An alternative way to show the effects of our improved approximation is to define 
an effective velocity dispersion, 

° 2 ef fA k ^,0) = k2{ ~ 2 _ 1)2 HGa(k;v,0)), (59) 
which, in the high-rc limit, reduces to 0^,(77,0) defined in Eq. (49). In figs. 12, 13 



we 



plot cTe//,i ( on t ne left) and c^j 2 (on the right) as a function of k: the line code is the 
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Figure 14. er^.. 1 (z) as a function of redshift for three different momentum scales. 
Line-codes as in Fig. [9] 



same as in figs. [9] and 10 
Finally, in figs. 



14 



and 



15 



values of the momentum: k = 



we plot a _ X £ „ CIS cl function of redshift for three fixed 



e//,o 

).5, 0.15, 0.001 h/Mpc. 
Our results show that, for z ~ 2, the subleading effects neglected in the Crocce- 
Scoccimarro approximation start to play a relevant role. For the density propagator 
G\ they are larger that 1% for k > O.lOh/Mpc at z = and for k > 0.25h/Mpc at 
z — 1. For the velocity propagator G 2 the effect is stronger. It is larger than 1% for 
k > 0.07h/Mpc at z = and for k > 0.16h/Mpc at z = 1. These effects should clearly 
be taken into account in a computation aiming to reproduce the BAO power spectrum 
at the percent level. At k — 0.2h/Mpc, that is, well inside the BAO range of scales, 
and at z — 0, the deviation from the Crocce-Scoccimarro resummation is 4.0% for G\ 
and 7.8% for G 2 . 

The sign of the correction is the same in both approximations considered for the 
non-linear PS {i.e. 1-loop and TRG): the propagators are less damped than the Crocce- 
Scoccimarro one. This can be easily understood analytically by looking at the difference 
between Eq. (38) and Eq. (49). In the former, the large-fc damping is modulated, via 
cr 2 , by the linear PS, P°, while, in the latter, it is modulated by the velocity-velocity 
component of the non-linear PS, i.e. by P^, see Eq. (50). Now, unlike the non-linear 
density-density PS, P-Jj, which is enhanced w.r.t. the linear one, receives negative 
corrections, and is therefore smaller than P° at intermediate and large fc's (see, for 
instance, [HI US])- Moreover in Eq. (49) the unequal-time cross-correlator appears, which 
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Figure 15. Same as Fig. 14 for a^ff 2 (z) 




Figure 16. Contributions at z — to non linear term of Eq. (58) for the density 
propagator (left panel) and velocity one (right panel). They are obtained by using 
1-loop approximation for the non linear PS. The solid continuos lines give us the whole 
/ E ac u c value, while the dotted green, the short-dashed magenta, the long-dashed cyan 
and the dash-dotted orange ones represent the P^ 1 , P^, P%1 and contributions 
respectively. 
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Figure 17. Same as Fig. 16 where now we consider the TRG approximation for the 
non linear PS. 



is further suppressed w.r.t. the equal-time one. As a consequence, we get cr^ < a 2 and 
therefore a smaller damping at large k. On the other hand, for small and intermediate 



momentum values, the non linear part of Eq. (58) depends also on Pf/, Pfj f° r the 
density propagator and on P%{ for the velocity propagator. In figs. 16 and 17 we 
plot the contributions to J ds S E 



j ac u c evaluated at redshift zero for 1-loop and TRG 
power spectra respectively. The solid continuos lines gives the whole / dsT, ac u c value, 
while the dotted green, the short-dashed magenta, the long-dashed cyan and the dash- 
dotted orange ones represent the Pq, P^ 1 , P^ 1 and P^ 1 contributions respectively. Notice 
that, also in this range of scales, the velocity- velocity component gives the dominant 



contribution to Eq. (58). In particular the non linear term that depends on the density- 



density component of the PS is positive but subdominant. Moreover, comparing the two 



figs. [16] and [17] it is clear that the unphysical behavior of the 1-loop PS approximation 
discussed above translates into an underestimation of the non linear contributions to 



Eq. (58). 



7. Discussion and conclusions 

In this paper we have extended the computation of the non-linear propagator pioneered 
in CS, which was based on the resummation of the chain- diagrams at all orders in PT. 
We have taken into account new contributions, obtained by replacing the linear PS 
appearing in the chain-diagrams by the non-linear PS. We have proved the remarkable 



property (see Appendix A), that this wider class of renormalized chain-diagrams can be 
exactly resummed in the large k limit. In the same spirit of CS, we required that PT is 
recovered in the k — > limit, which implies taking into account diagrams not belonging 
to the renormalized-chain diagrams class. 

The resummation of this extended class of diagrams is greatly simplified by the use 
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of the time evolution equation for the full propagator, Eq. (58). Instead of dealing with 
a complex diagrammatic analysis, the task is reduced to the solution of a differential 
equation for the full propagator. The crucial element in this equation is the self-energy 
E a fe. Approximating it with the 1-loop self-energy gives the CS result. On the other 
hand, the renormalized chain-diagrams resummation is achieved by considering a E a & 
which is still formally 1-loop, but with the linear PS replaced by the non-linear one, see 
Fig. |4| We have tested two different approximations for the non-linear PS: the 1-loop 
approximation, and the result of the TRG evolution [9]. The former gives results that 
depend sensibly on the UV cutoff in the loop integral. This is to be expected since, at 
low redshift, the 1-loop PS becomes unreliable, and even negative, al large momenta. 
On the other hand, the PS from the TRG evolution does not suffer from UV problems. 
Moreover, as discussed in detail in [5], the solution of the TRG equations for the PS is 
formally a 1-loop expression, in which the linear PS' are replaced by non-linear ones, 
and is therefore fully consistent with the spirit of our treatment for the propagator in 
this paper. 

The numerical results show that the new effects are quite relevant in the BAO scales, 
where they are in the few percent range at z=0. They should therefore be taken into 
account in computations of the PS in the BAO range based on the use of renormalized 
propagators, such as RPT. Indeed, in [7J, the effect of next-to leading order corrections 
to the propagator was advocated in order to reconcile the RPT results on the PS with 
N-Body simulations. The authors correctly identified the effect of these corrections 
with a renormalization of the linear PS which, through the quantity a 2 , modulates the 
Gaussian decay of the propagator at large momentum. However, instead of performing 
an explicit computation as the one presented in this work, they implemented an ad 
hoc procedure, by replacing the linear PS with the non-linear one as obtained in the 
halo model [26, 27J. This is inconsistent, since, by doing so, one is using the density 



PS, whereas the large-/c limit resummation involves the velocity PS, see Eq. (48). As 
a result, the procedure illustrated in [7J leads to a wrong prediction on the sign of the 
corrections to the CS propagator induced by the subleading corrections. 

We stress that our conclusion that the effect of these corrections is to enhance the 
propagator w.r.t. the CS result is by no means based on a particular approximation 
of the non-linear PS. Different choices for the latter may give different results on the 
size of these corrections, but the sign is only determined by the assumption that the 
non-linear velocity PS is smaller than the linear one, which is verified in all consistent 
approximations (see, for instance, [9| [13] ) 

A careful reconsideration of the comparison with N-body simulation is therefore 
needed, both for the density and for the velocity propagators. Should the discrepancy 
mentioned in [21 [7] persist, it would imply that other effects should be taken into account. 
One possibility would be to include diagrams not belonging to the renormalized chain 
class. At 2-loop order, it would mean to include diagrams VII, VIII, and IX in Fig. [7j 
also at large k. Notice that these contributions do not exponentiate in the large-/c limit 
or, equivalently, they break the factorization property of Eq. (34) and, consequently, 
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they give a propagator which deviates from the gaussian decay form found by CS at 
large k. Another possible reason for the discrepancy could be the effect of small scale 
non-linearities, which translates in a non-vanishing velocity dispersion at intermediate 
scales see, for instance, [28]. The inclusion of this effect in the computation of the 
(resummed) propagator will be analyzed elsewhere. 

The approach followed in this paper is based on the use of the exact evolution 
equation (30). On the same spirit, an evolution equation for the PS can be written, 
and its results can be compared with alternative approaches presented in the recent 
literature [9j [121 02] • This will be the subject of a forthcoming publication. 
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Appendix A. Factorization at large momentum 



In this appendix, we will show that, at large momentum, Eq. (34) holds, that is 
£ rdsZ£-i\k; Va ,s)G < £(k; s, Vb ) 

~'^0 lb 



; Va,S)U c . 



(a.i; 



The first point to notice is that, as shown by CS in p], the leading contributions at 
large k and at a fixed loop order n goes as k 2n , and is given by the chain- diagrams of 
Fig. [3| These diagrams are such that the 2n propagators are lined up in a single chain 
carrying the momentum k, and any of the n power spectra, carrying a lower momentum 
qi, is connected to the propagator chain by both its legs. Each of the 2n vertices now 
contributes a factor 



,, , x largefc 1 k • q 

w c 7ad>(k, -q,q - k) — Y -^-d ab , 

Z q 



(A.2) 



where the u c comes from the PS (see Eq. (A. 3) below), giving the above mentioned 
0(k 2n ) behavior. Since the linear propagators are momentum- in dependent, the 
integrals over the loop momenta decouple one another, each one giving a contribution 
proportional to the 1-loop 'self-energy', 

laiddAK qi, -k - (li)g dihi (s ai - SbJ'Jh^ik + q t , -q i; -k) 



large k 



_ k 2 a 2 e s ai +s b . Qaibi ( Sai _ s ) 



where 



a 



q 2 



(A.3) 



(A.4) 



and we have used Eq. (A.2). 



Thanks to the composition property of the linear propagators 

gac(Va - Vc)9cb(Vc - Vb) = 9ab(Va ~ Vb) , 



(A.5) 



the chain of linear propagators emerging in the large k limit of (A.3) combine into a 
single one, g a b(Va — Vb), independent of the intermediate times. 

In order to discuss the time integrals, we consider a generic n-loop contribution 



to the sum (A.l) (see Fig. Al for a 5-loop example). The 'self-energy' is a n — j loop 
quantity and the propagator a j-loop one. Fixing the intermediate time s, the 'self- 
energy' diagram has thus 2(n — j) — 2 intermediate times, while the propagator has 2j. 
The 'self-energy' time integrals give 

a, :■ ~r .' 

&2(n-i)-2 



dti / dt 2 ■ ■ ■ / dfo fn -n-* e^+^-i 

J S J S 



(2n-2j-2)\ 



(A.6) 
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Figure Al. A 5-loop contribution to the sum of Eq. ( A.l I 



On the other hand, the propagator time integrals are 



dn dr 2 - ■ ■ I dr 2j e^k=i 

JVb J^b ->Vb 



(e s - e Vb ) 2j 

m 



(A.7) 



Multiplying (A. 6) by (A.7) and by the remainig e s time factor, and then integrating 

(A.8) 



over s from r\b to r) a , as in the LHS of (A.l) one gets the time coefficient 



(2n-l)! 



which is independent of j, i.e. is the same for any term in the sum in Eq. (A.l) 



and depends only on the total loop order, n. The sum over j ensures that the n power 
spectra are attached to the propagator chain in all possible ways, i. e. that all the pairings 
between the 2n vertices are taken into account. There are (2n — 1)!! such pairings, so, 



using Eqs. (A.3) and (A.8), we obtain the LHS of (A.l) in the large momentum limit 



(2n - 1)! 
1 



(2n-l)\\(-k 2 a 2 ) n 



(n-1)! 
1 

(n - 1)\ 



2 2 ' 



o'la 



oVb\2 



-k 2 a 2 



2 

, < <ia - e 



n— 1 



Va _ p^O 2 " 1 "" 1 



{-k 2 a 2 )e^[e^ - e % ) 
ds T^(k\ r) a ,s)u c . 



Va 
Vb 



(A.9) 



The contribution to the propagator from chain-diagrams at (n — l)-loop order can be 
computed straightforwardly, by using the same properties considered above pQ. The 
integration over the 2(n — 1) intermediate times gives 



Va rti 

dt\ / dt 2 

Vb JVb 



Vb (2n-2)! 



(A.10) 



Since there are (2n — 3)!! chain-diagrams at n — 1 order, the propagator in the large k 
limit reads 



large k 



1 



(n-1)! 



2 2 



o'la 



,Vb\2 



-i n-1 



-KG 



g a b(va-Vb), (A.ll) 



which, comparing with the last line of ( |A.9 ), proves Eq. (A.l). 



The factorization property, Eq. (A.l), is a general property of chain-diagrams, 



which does not rely on the fact that the PS considered above is the linear one: 
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P a b(q; s a , s b ) = P°(q)u a u b . Indeed, in the remaining part of this Appendix, we will 
show that the factorization still holds if one considers a generic form for the PS, i.e. 

Pab{q; s a , s b ) = Pba(q; s b , s a ) , (A.12) 

and therefore it holds if one renormalizes the linear PS by including non-linearities in 
different consistent approximations, such as, perturbation theory, TRG, and so on. 

The starting point is to realize that the LHS of Eq. ( A.l ) is obtained by pairing in 
all possible ways the 2n vertices - including the extremal one at time rj a - connected 
by the chain of 2n propagators. Moreover, the large momentum property of the vertex, 
Eq. (A. 2), still holds of it is contracted by a generic vector A a , giving 

lk- q 



A c j acb (k , -q , q - k) 1 ^ fc A 2 



Therefore, in the large k limit, we have 



2 q 2 



Sab ■ 



(A.13) 



n-1 



'la 



j=Q Jr >b 



large k 



Va 



'lb 



Si 



dsi / ds2 



>ib 



S2n-2 



ds 



2n-l 



'lb 



n« =1 A* ) x 



2n pairings 



-f22(9i; Va, S ai ) P22(q n ] Sa 2n -2, S a 2 „-i, 



Ql 



q 2 
in 



(A.14) 

Notice that the time-integrand is, by construction, invariant under the exchange of any 
of the 2n — 1 variables, Sj -H- Sj, therefore we will use the property 



Va 



in, 



si 



dsi / ds2 



'ib 



sjv-i 



ds N jF[si, ■ ■ ■ ,s N ] 



'ib 
1 

M 



Va 



'lb 



dsi ds 2 



in, 



Va 



dsjv^fsi, • • • ,s N ] , (A. 15) 



>ib 



where the function J-[sx, • • • , sn] is totally symmetric. Eq. (A.14) can then be rewritten 

as 



'-k 2 



(2n - 1)! 



I la 



ds\ I ds 



2 • • • 



vi< 



in. 



ds 2n -i [ I njUAi ) x 



2n pairings 



f22(gi; Va, S ai ) P22{qn, S a2n _ 2 , S a2n _A 



Qi 



-k' 2 



'la 



ds e Va 



in, 



d 3 q 



P22{q;Va,s) 



X 



(2n - 2)! 



'la 



'ib 



'la 



ds<? 



'ib 



'la 



ds 



2n~2 



Vb 



Ai ) x 



2n—2 pairings 



^22(9i; 



Pllifln-l, s a2n-z, S ain-2, 
n 2 

Hn-1 



(A.16) 



Now the upper indices count the number of non-linear PS in a given contribution to £ and G. 
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which, multiplied by g a b(Va — T] b ), gives the RHS of (A.l) (see also (A. 3) and (A. 11)). 



Appendix B. Time-Remornalization Group equations for the non linear PS 

In this Appendix we briefly review the Time-Remornalization Group approach (TRG) 
introduced in [5] to compute the non linear PS. Applying the equation of motion in 
Eq. ^ to the (equal-time) PS, the bispectrum, 

(<Ax(k, ??Vb(p, v)<Pc(SL, V)) = M k + P + q) B abc (k, p, q; rj) , (B.l) 

and to the higher order correlators, one gets an infinite system of coupled differential 
equations. Truncating the hierarchy by setting the trispectum (i.e. the connected four- 
point function) to zero, one is left with the closed system 

d v P ab (k ; 77) = -Q ac P cb (k ; 77) - Q bc P ac (k ; 77) 

+ e v J d 3 q [7 ac d(k, -q, q - k) B bcd (k, -q, q - k; 77) 

+B acd (k, -q, q - k; 77) 7 fecd (k, -q, q - k)] , 



d v B abc (k, -q, q - k; 77) = -Q ad B dbc (k, -q, q - k; 77) 

- tt bd B adc (k, -q, q - k; 77) 

- Q cd B abd (k, -q, q - k; 77) 

+ 2e' 7 [7ad e (k, -q, q - k)P db (q ; r])P ec (k - q ; 77) 
+ 76<fe(-q, q - k, k)P dc (k - q ; 77)P ea (k ; 77) 
+ 7cde (q-k, k, -c l )P da {k;r ] )P eb {q;r ] )] . (B.2) 



The formal solution of the system (B.2) is given by 

P ab ik ; rj) = g ac (r], 0) ^(77, 0)P cd (k ; 77 = 0) 



+ / drj'e* 1 J d q g a e(v,v')9bf(v,v') 

x hecd(k, -q, q - k) B fcd (k, -q, q - k; 77') 

+ 7/cd(k, -q, q - k) B ecd (k } -q, q - k; 77')] 



B abc (k, -q, q - k; 77) = 

9ad(v, 0)5'b e (?7, 0)g c f(v, 0)B def (k, -q, q - k; 77 = 0) 

rv , 

+ 2 / dr]'e v g ad (r) 7 r))g be (r) 7 rj')g cf (r], 77') 
J 

x [ld 9 h(k, -q, q - k)P e9 (q ; 7]')P fh (q - k ; 7]') 
+ 7egh(-q, q - k, k)P /9 (q - k ; r]')P dh (k ; 77') 
+7 /j/ft (q-k, k, -q)P dg {k- v ')P eh (q ]V ')} , (B.3) 

which shows that the non-linear PS in this approach is given by a formally 1-loop 
expression in which non-linear PS' replace the linear ones. In this respect, the TRG 
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approach is fully consistent with the computation of the propagator presented in this 
paper, which is based, as well, on improving over the CS approximation by replacing 
linear PS' with non-linear ones. 
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